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Abstract 

It is shown that the infall of collisionless dark matter onto isolated galaxies 

produces a series of caustic rings in the halo dark matter distribution. The 

properties of these caustics are investigated. The density profile of the caustic 

is derived for a specific case. Bumps in the rotation curve of NGC 3198 are 

interpreted as due to caustic rings of dark matter. 
PACS numbers: 95.35.+d, 98.35.Gi 
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There are compelling reasons to believe that the dominant component of the dark matter 
of the universe is non-baryonic collisionless particles |1[ . The leading candidates are axions, 
WIMPs and massive neutrinos. The word "collisionless" signifies that the particles are so 
weakly interacting that they have moved purely under the influence of gravity since their 
decoupling at a very early time (of order 10~ 4 sec for axions, of order 1 sec for neutrinos and 
WIMPs). In the limit where the primordial velocity dispersion of the particles is neglected, 
they all lie on the same 3-dim. 'sheet' in 6-dim. phase-space. Their phase-space evolution 
must obey Liouville's theorem. This implies that the 3-dim. sheet cannot tear and hence 
that it satisfies certain topological constraints. 

Let us assume that collisionless dark matter (CDM) exists. Usually, CDM means 'cold 
dark matter' (e.g., axions or WIMPs) but because massive neutrinos (m„ >> 10 _1 eV) 
behave, from the point of view of this paper, in a similar way we include them in our 
definition as well. Because their phase-space sheet cannot tear, CDM particles must be 
present everywhere in space, including specifically intergalactic space. The space density 
may be reduced by stretching of the phase-space sheet but it cannot vanish. Moreover, the 
average space density is recovered as soon as the average is taken over distances larger than 
the distance CDM may have locally moved away from perfect Hubble flow. In a region 
which is sparsely populated with galaxies, this distance is much smaller than the distance 
between galaxies. The implication is that isolated galaxies are surrounded by unseen CDM 
and hence, because of gravity, CDM keeps falling onto isolated galaxies continuously from 
all directions. If the galaxy joins other galaxies to form a cluster, infall onto the galaxy gets 
shut off because of lack of material but infall onto the cluster continues assuming that the 
cluster is itself isolated. In an open universe (fi < 1), the infall process eventually turns off 
because the universe becomes very dilute. However, even if our universe is open, we are far 
from having reached the turn-off time. 

It has been shown that, under a wide range of circumstances, CDM infall onto an 
isolated galaxy produces a halo whose density falls off like \ for large r, where r is the 
distance to the galactic center. Such a halo implies a flat galactic rotation curve 0. It has 
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also been shown 0] that the angular momentum carried by the CDM particles has the effect 
of depleting the inner halo and of making the halo contribution to the rotation velocity 
vanish at r — > 0, thereby introducing an effective core radius. Both the approximately flat 
rotation curves and the presence of effective core radii are consistent with observation. CDM 
infall causes peaks in the velocity distribution of dark matter particles on Earth The 
present paper focuses on the related phenomenon of ring caustics in the halo density. The 
existence of these caustics was noted in ref. (§] but their properties were left unexplored. 

Consider the infall of CDM onto an isolated galaxy. Let's first neglect the velocity 
dispersion of the infalling particles. In practice it is sufficient that the velocity dispersion 
of the infalling matter is much smaller than the rotation velocity of the galaxy but it is 
convenient to consider the extreme case of zero velocity dispersion first. Consider the time 
evolution of all CDM particles which are about to fall onto the galaxy for the first time in 
their history at time t. For the sake of definiteness, we may consider all particles which have 
zero radial velocity (f = 0) for the first time then. Such particles are said to be at their 'first 
turnaround'; they were receding from the galaxy as part of the general Hubble flow before 
t and will be falling onto the galaxy just after t. They form a closed surface, enclosing the 
galaxy, hereafter called the turnaround 'sphere' at time t. (The present turnaround sphere 
of the Milky Way galaxy has a radius of order 2 Mpc.) The turnaround sphere at time t falls 
through the central parts of the galaxy at a time of order 2t. Particles falling through the 
galactic disk (assuming the galaxy is a spiral) get scattered by an angle A6 ~ 10^ 3 by the 
gravitational fields of various inhomogeneities such as molecular clouds, globular clusters 
and stars |J. However most particles carry too much angular momentum to reach the 
luminous parts of the galaxy and are scattered much less. Because the scattering is small, 
the particles on the turnaround sphere at time t, after falling through the galaxy, form a 
new sphere which reaches its maximum radius R' at some time il . The radius R' at the 
second turnaround is smaller than the radius R at the first turnaround because the galaxy 
has grown by infall in the meantime. The sphere continues oscillating in this way although 
it gets progressively fuzzier because of scattering off inhomogeneities in the galaxy. 
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From a topological viewpoint, each time the sphere falls through the galaxy it turns itself 
inside out and hence there is a ring of special space-time points associated with each fall 
through the galaxy. This ring may be defined as the locus of points which are inside the 
sphere last, just before the sphere completes the process of turning itself inside out. Fig. 
I illustrates this. It shows successive 2-dim. cross-sections of a sphere falling in and out. 
After the sphere has crossed itself, as in frame (c), there is a donut-shaped volume of points 
which are still inside the original sphere. The ring is where this donut-shaped volume ends in 
space-time. Because CDM falls in continuously, the ring in space-time associated with one 
infalling sphere is in fact a persistent feature in space. For an arbitrary angular momentum 
distribution on the turnaround sphere, the ring is a closed loop of arbitrary shape. However, 
if the angular momentum distribution is dominated by a smooth component which carries 
net angular momentum, the ring resembles a circle. If there is no angular momentum at 
all, the ring reduces to a point at the galactic center. As we will soon see, the ring is the 
location of a caustic with strong density contrast. There is one caustic ring due to CDM 
particles falling through the galaxy for the first time, a caustic ring of smaller radius due 
to particles falling through for the second time, a yet smaller ring due to particles falling 
through for the third time, and so on. The caustic rings move slowly in space, generally in 
an outward direction, as the properties of the infalling CDM particles (in particular, their 
turnaround radius and angular momentum distribution) change. 

Let's obtain the density profile of the caustic in a particular case of axial symmetry 
about the z-axis and parity (z— >-z). The symmetry assumptions will force the caustic to lie 
on a circle in the z=0 plane. Consider the time evolution of CDM particles initially located 
with uniform density on a turn-around sphere of radius R and with the angular momentum 
per unit mass distribution £(f) = ufx [z x r) where r — Rf. Thus the turnaround sphere 
is assumed to be initially rotating about the z-axis with angular velocity u as if it were 
a rigid body. Let's assume further that the particles fall into the spherically symmetric 
potential: V(r) = v^ ot \n(R/r), which is such that it produces a flat rotation curve with 
rotation velocity v rot . One readily finds that a continuous flow of particles falling through 
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under the stated conditions produces a density distribution which has a caustic ring in the 
z = plane of radius a given by: 
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a v rot p \n(R/a) = wR 2 (l + 0(— )) . (1) 

The caustic is the locus of points where the particles with the largest amount of angular 
momentum (£ m ax = ^R 2 ) are at their distance of closest approach to the galactic center. A 
calculation shows that, to leading order when a « R, the density distribution is given by: 

dM 2 1 
d{a;p,z) = —— . (2) 

dfldt V ^/( r 2 _ fl 2)2 + 4a 2 z 2 

where (p, z, 6) are cylindrical coordinates, r = \/ p 2 + z 2 , JfL is the rate at which mass falls 



in per unit time and unit solid angle, and v = v rot y2\n(R/a) is the velocity of particles near 
r = a. From Eq. (2), one readily finds the behaviour of the density near the caustic: 

dM 1 

d{a;p,z)~ — — 3 

aildt v a a 



where a = \J (p — a) 2 + z 2 is the distance to the ring. Of course the caustic singularity 
is only perfectly sharp in the limit where the velocity dispersion of the infalling CDM is 
zero. If the velocity dispersion is Av, the caustic singularity spreads over a distance of order 
Aa = R^. 

For an isolated galaxy, there is a ring of radius a\ due to particles falling through for 
the first time, a ring of radius 02 due to particles falling through for the second time, and 
so on. Thus the quantities R, a, v, d and J|4 carry an index n = 1,2,3.... One caustic 
ring is associated with each pair of velocity peaks [||,U in the halo distribution. Eq. (2) 
implies j^^- = r 2 d n (0;r). This relates the value of the prefactor ^2iv~ ^ or eacn m f a ^ 
to the density d n (0;r) the associated pair of velocity peaks contributes in the limit of zero 
angular momentum (a = 0). Estimates of d n (0,r) can be found in ref. Q for the case of 
self-similar infall. 

The infall is called self-similar if it is time-independent after all distances are rescaled 
by the turn-around radius R(t) at time t and all masses are rescaled by the mass M(t) 



interior to R(t). In the case of zero angular momentum and spherical symmetry, the infall 
is self-similar if the initial overdensity profile has the form = (If 1 ) 6 where M and e are 
parameters 0. e must be in the range < e < 1. The rotation curve is flat if0<e<2/3 
0. The infall model was generalized in ref. to include the effect of angular momentum. It 
was found that self-similarity requires the angular momentum distribution E(t) to have the 
time-dependence £(t) = j — ^ , where j is a dimensionless and time-independent distribution. 
A good agreement of the self-similar model with the properties of our own galaxy was found 
0] for parameter values of order e = 0.25, j = 0.2 and h = 0.7 where j is the average of 
the j distribution and h is the Hubble rate in units of 100 km/sec. Mpc . Using the model 
of ref. m, the following values for the radii of the first five rings are obtained: 

K : n = 1, 2...S} * (39, 19.5, 13, 10, 8)kpc (^) (^) ( J^) (4) 

where j max is the maximum j value of the j distribution, and where the value e = 0.3 was 
used. For e = 0.2, a\ ~ 3Qkpc {^f^j {jjf) ( 220km /s ) ^ e ra ti° s a n/^i are almost the 
same as in the e = 0.3 case. 

Luminous rings surrounding galaxies have been observed and such rings may be related 
to the caustics described here. However, luminous matter is presumably baryonic and the 
behaviour of baryons is more complicated than that of CDM particles. Infalling baryons 
do not necessarily behave in a collisionless manner and they may easily get stripped off 
the CDM flow. In contrast, the conclusion that isolated galaxies are surrounded by caustic 
rings of dark matter appears unavoidable if CDM exists. Because of this, we will limit 
ourselves in this paper to the observational implications of CDM rings which follow from 
their gravitational fields. CDM rings may nonetheless have a baryonic component. These 
baryons may have the same phase-space distribution as the CDM or they may have accreted 
onto the ring from neighboring space. 

The Newtonian gravitational force per unit mass due to the density distribution of Eq. 
(3) is 

F n (a) = (5 ) 

On 
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for small er, where C n = ^jt~v~ = r ' 2 ^n(0; r). Hence there is a discontinuity in the rotation 
velocity 

A n Vrot _ 1 d n (0;a n ) 
v rot 2 d(a n ) 

directly across the n th caustic, in the ideal case where this caustic is perfectly sharp and lies 
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in the galactic plane. In Eq. (6), d(r) = 4 ^, 2 is the total density at r in the limit of spherical 
symmetry and perfectly flat rotation curves. The ratios /„ = \ dn }^ were f° un d to be of 
order (13, 5.5, 3.5, 2.5, 2, ...)10~ 2 in the self-similar infall model with e = 0.2. (For r << R 
and small n, the f n are nearly r- independent.) There are of course a number of effects that 
will smooth out sudden variations in the measured rotation velocities. One effect is that the 
CDM ring may be some distance away from the galactic plane where the rotation velocities 
are measured. Secondly, the measured rotation velocities are spatial averages over some 
distance across the galactic plane. Thirdly, the CDM ring may be fuzzy because of velocity 
dispersion as mentioned earlier. Thus the discontinuities of the ideal case are smoothed out 
into bumps. The bumps occur in the measured rotation curve near the location of CDM 
rings if the latter happen to be close to the galactic plane. 

Galactic rotation curves often do have bumps. Of special interest here are those bumps 
which occur at radii larger than the disc radius because they cannot readily be attributed to 
inhomogeneities in the luminous matter distribution. Consider the rotation curve of NGC 
3198 0, one of the best measured and often cited as providing compelling evidence for 
the existence of dark halos. It appears to have bumps near 28, 13.5 and 9 kpc, assuming 
h = 0.75. Although the statistical significance of these bumps is not great, let's assume for 
the moment that they are real effects. Note then that their existence is inconsistent with 
the assumption that the dark halo is a perfect isothermal sphere. On the other hand, the 
radii at which the bumps occur are in close agreement with the ratios: 

{a n /a n _i : n = 2,3,4,5} = (0.49, 0.67, 0.76, 0.81) (7) 

predicted by Eq. (4) assuming that the bumps are caused by the gravitational fields of the 
first three caustic rings of NGC 3198. Since v rot = 150 km/sec, we find that j max = 0.28 in 
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this case if e = 0.3. The uncertainty in h drops out. A fit of the infall model to our own 
galaxy [|[ produced j ~ 0.2 for e = 0.2 to 0.3. If the turnaround sphere is taken to be rigidly 
rotating, one has j max = -j. Thus the values of j for our own halo and that of NGC 3198 
are found to be similar. 

The ratios of caustic ring radii given in Eq.(7) are characteristic of the t-dependence in 
the angular momentum distribution l(f, t) = j(f)R(t) 2 jt . The main reason for using this 
ansatz in ref. 0] was that it produces exact self-similarity in the time evolution of the halo. 
However, there is a broader justification for a time- dependence close to the one of the ansatz. 
Consider particles which at some early initial time t{ are at a distance from the center 
of a large overdensity which will grow into a galactic halo. These particles have magnitude 
of angular momentum with respect to the overdensity's center: /(r*j) = r^Vix^i) where Vi± 
is the magnitude of the component of the initial velocity Viifj) transverse to r*j. Because 
v(0) = 0, it is reasonable to assume Vi±(fi) is proportional to when comparing values of 
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Ti which are of the same order of magnitude. In that case, Z(r*j) ~ r\ ~ M? ~ i?(rj) 2 /t(rj) 3 
where Mj is the mass interior at time tj to the sphere of radius rj, R(ri) and t(rj) are the 
turn-around radius and turn-around time of particles initially at radius r^, and where we 
used GMi = n 2 R(ri) 3 /8t(ri) 2 . Except for the relatively slowly-varying factor of l/ts, this 
is the time dependence which leads to self-similarity. At any rate, since the gravitational 
field around a large isolated overdensity is approximately spherically symmetric and hence 
angular momentum about the center of the overdensity is approximately conserved, the 
study of bumps in rotation curves may inform us about the peculiar velocities associated 
with primordial density perturbations and thus constrain theories on the origin of these 
perturbations. 

Finally, let's discuss the implications for our own halo of the assumption that the caustic 
rings lie in the galactic plane. As an example, we will use the model parameters e = 0.25, 
jmax = 0.265 and h = 0.7. In this case the Milky Way caustic rings are at the radii: 41, 20, 
13, 10, 8.0, 6.7, 5.8, 5.1, 4.5, 4.0, 3.7, 3.4, 3.1 ... kpc. We at 8.5 kpc would be between the 
4th and 5th ring. There is evidence for the 6th through 13th caustic rings in that there are 
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sudden rises in the inner rotation curve published in ref. || at radii very near those listed 
above. Similarly there is some evidence for the 2d and 3d caustic rings in the averaged 
outer rotation curve published in ref. ||. This will be discussed in detail in a future paper 
|I(J . The 4th ring would have passed by us approximately 900 million years ago assuming 
that our distance to the galactic center is and has always been 8.5 kpc. The 5th ring would 
pass by us approximately 350 million years from now. It is conceivable that the passage of 
caustic rings caused cataclysmic events in our past. If the angular momentum distribution 
is the one ^ = j6>( -^ a3: ~ J L characterizing a turn-around sphere which is initially rigidly 
rotating, the contributions of the first eight incoming peaks to the local halo density would 
be approximately (0.4, 0.9, 1.7, 5, 13, 3, 2, 1.4) lO -26 ^- whereas the averages || over all 
locations at the same distance from the galactic center as us are approximately (0.4, 0.8, 1.2, 
2, 2.7, 1.6, 1.2, 1.0) 10~ 26 ^3. The 4th and 5th velocity peaks are considerably enhanced 
because of our proximity to the corresponding rings. In the model, the two flows associated 
with the 4th ring have velocity vectors (2100 ± 260z)km/s in the Sun's rest frame whereas 
those associated with the 5th ring have velocity vectors (2300 ± 160f)km/s. Here, z is the 
direction perpendicular to the galactic plane, f is the radial direction in the galactic plane 
and is the direction of rotation in the plane. If the caustic rings lie in the plane, the 
local halo density is boosted from approximately 0.5 lO" 24 ^-, which is the usual estimate 
in the case of spherical symmetry, to approximately 0.9 lO -24 ^-. These results are relevant 
to the axion [|ll]] and WIMP [|l!| dark matter searches which in fact provided the original 



impetus for this work. It is conceivable that these experiments will measure some day the 
contributions of the velocity peaks to the local density and thereby provide us with detailed 
information about the structure of our halo. 
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FIGURES 



a) b) 




FIG. 1. Infall of all particles on a turnaround sphere. The sphere has net angular momentum 
about the vertical axis. It crosses itself between frames b) and c). The cusps in frames d) and e) 
are at the intersection of the ring caustic with the plane of the figure. 
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